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Abstract 

We examine spectra of Dirac operators on compact hyperbolic surfaces. Particular 
attention is devoted to symmetry considerations, leading to non-trivial multiplicities 
of eigenvalues. The relation to spectra of Maafi-Laplace operators is also exploited. 
Our main result is a Selberg trace formula for Dirac operators on hyperbolic surfaces. 
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1 Introduction 



Trace formulae play a prominent role in spectral geometry and in quantum chaos. In spec- 
tral geometry they relate spectra of certain (pseudo-) differential operators on manifolds 
to the geometry of that manifold. The prime example of such a trace formula, the classical 
Selberg trace formula |Sel56j , is concerned with the Laplace-Beltrami operator on a hyper- 
bolic manifold, i.e., a Riemannian manifold with constant negative sectional curvatures. 
This trace formula establishes a connection between the spectrum of the Laplacian and 
the length spectrum of closed geodesies on the manifold. 

In a semiclassical context spectra of semiclassical operators are related to the periodic 
orbits of an associated Hamiltonian flow. In extensions of this procedure to operators 
acting on sections in vector bundles over the relevant phase space, it turned out that 
the Hamiltonian flow does not provide the entire input on the classical side of the trace 
formula (see, e.g., |BK99j). In particular, for the Dirac operator on M 3 the corresponding 
classical dynamics are given in terms of skew-product flows over Hamiltonian flows, with 
spin precessions as their cocycles. A similar result is found in the case of a trace formula 
for the Dirac operator on a graph BH03,. Very recently, in the context of quantum 
ergodicity Jakobson and Strohmaier [,IS06| observed that (the square of) a Dirac operator 
on a compact Riemannian manifold is related to the associated frame flow, which is a 
geometric analogue of the skew-product flows in the semiclassical setting. 

Here we introduce Dirac operators on compact hyperbolic surfaces, investigate their 
spectra in terms of symmetries, and finally develop the related Selberg trace formula in 
a classical setting. This is a direct approach employing Green's functions and point- 
pair invariants, and therefore complements previous work B095J, where trace formulae 
for positive square-roots of squared Dirac operators were determined using Lie-algebraic 
methods. On the way we exploit the well-known relation of the squared Dirac operator to 
a Laplacian that contains a coupling to a magnetic field (see, e.g., |Pnu94j ) . In particular, 
we show that spectra of these MaaB-Laplacians, when they correspond to odd weights, 
possess multiplicities of at least two. This is a consequence of Kramers' degeneracy in the 
spectrum of the Dirac operator. The geometric side of the trace formula for the Dirac 
operator is primarily determined by the closed geodesies on the surface. In addition, we 
identify traces of the frame flow and its lift into the spin structure, which provide suitable 
phase factors associated with the closed geodesies. 

The paper is organised as follows: In Section |2] we review the definition of a Dirac 
operator on a Riemannian manifold and explicitly carry out the constructions in the case 
of the hyperbolic plane. Section Elthen is devoted to an identification of spinor bundles over 
compact hyperbolic surfaces via factors of automorphy. We then investigate symmetries 
of the Dirac operator on a compact surface and study their influence on the spectrum in 
Section |U After having introduced point-pair invariants in Section El we calculate traces of 
Hilbert- Schmidt operators to arrive at the desired trace formula in Sectional Finally, in 
Section^ we discuss some applications. 
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2 The Dirac operator 



A Dirac operator is a first order, elliptic differential operator acting on sections in a spinor 
bundle over a Riemannian manifold. For its construction one first needs to supply an 
orientable, n-dimensional Riemannian manifold (M, g) with a spin structure. To this end, 
let FM be the SO(n)-principal bundle of oriented orthonormal frames of TM. A Spin(n)- 
principal bundle Q is called spin structure on M if there exists a principal-bundle morphism 
(p : Q —* FM that is equivariant with respect to the two-fold covering A : Spin(n) — > SO(n) 
(see e.g. |Fri97|lBau81| ). Moreover, in even dimensions n the complexified Clifford algebra 
<L c n is isomorphic to M(2 n / 2 , C), i.e., the algebra of complex 2 n//2 x 2 n//2 matrices. The action 
of M(2 n//2 , C) on C 2 ™ 2 defines the Clifford module A n . As Spin(n) is contained in A n 
is also a module for the spin group. This yields the spin representation p : Spin(n) — > 
Aut(A n ), which can be used to associate the spinor bundle S := Q x p A n to a given spin 
structure. 

The definition of the Dirac operator D rests on the canonical connection on Q that 
is given by the lift of the Levi-Civita connection on FM. Denoting the space of smooth 
sections of S by C°°(S), this connection induces a covariant derivative V s : C°°(S) — ► 
C°°(S®T*M). One then defines D = poio V s , where i denotes the canonical isomorphism 
between T*M and TM induced by the metric g, and \x is the Clifford multiplication of a 
vector field and a smooth section in S. 

The space C£°(S) of smooth, compactly supported sections in S can be turned into a 
pre-Hilbert space by introducing the inner product 



where d/x(m) denotes the volume form on M. The resulting Hilbert space will be denoted 
by L 2 (S'). The Dirac operator D is elliptic and essentially self-adjoint on C£°(S), thus its 
spectrum is real. Moreover, on a compact manifold M the spectrum is discrete. 

In the sequel we will focus on compact surfaces of constant negative curvature. Their 
universal covering space is the upper half-plane H 2 := {(x, y) £ IK. 2 1 2/ > 0} endowed with 
the Poincare metric g = y~ 2 ( dx <g> dx + dy <E> dy) of constant negative Gaussian curvature 
K = —1. The group of orientation-preserving isometries PSL(2,M) = SL(2, R)/{±Id} acts 
on H 2 via fractional linear transformations, 



where 7 = ( a c b d ) £ SL(2, M) is a representative for an element in PSL(2, M) and z = x + iy is 
viewed as a point in M 2 . Any compact hyperbolic surface can now be represented as T\EI 2 , 
where T C PSL(2,M) is a discrete, strictly hyperbolic subgroup (a cocompact Fuchsian 
group of the first kind). 

On the simply connected upper half-plane the relevant bundles are trivial, i.e., FM 2 ~ 
H 2 x SO(ra) and Q ~ H 2 x Spin(n). This allows an immediate construction of the Dirac 
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operator, which is explicitly given as the matrix-valued operator 

/ iyO + 8_1 

D = i . a o y d* y dy 2 

see, e.g., also [Pnu94j. The operators appearing in the off-diagonals have been introduced 
by MaaB and have been further investigated by Roelcke ( |Maa53t IRoe66t [Roe67j ) . In fact, 
they considered the operators 

d d k 
K fc := iy— + y— + - , 
ox oy 2 

. . d d k 

A fe := iy-z- - Vtt- + 77 • 

ox oy 2 

It is therefore natural to introduce the weighted Dirac operators 

D^i(_° At Y 2 ) • (") 

which reduce to the Dirac operator for k — 1. In physical terms the additional parameter 
k, which from now on will be called a weight, corresponds to a constant magnetic field on 
the surface. In normalised units the field strength is given by 

k- 1 / % 

B = . 2.2 

2e V ; 

On a compact surface the Dirac quantisation condition for the magnetic flux implies that 
k must be an integer. Therefore, from now on we will only consider fceZ. 

The weighted Dirac operator D& is closely related to the weighted MaaB-Laplacian, 

A k = -K k . 2 A k - |(1 - |) = y\£ 2 + *,) - iky£ , 

since 

Dk ~{ -A fc _ 2 -|(l-f) J • (2 ' 3) 

We remark that both the weighted Laplacians and the weighted Dirac operators are elliptic 
operators as their principal symbols do not depend on k. 

As yet we have neither specified spaces of sections in spinor bundles S, nor have we 
introduced Dirac operators on compact surfaces. It turns out that both purposes can be 
conveniently achieved by identifying suitable automorphic forms. We devote the following 
section to this task. 

3 Automorphic forms 

All relevant bundles over a compact surface r\H 2 correspond to trivial bundles over HI 2 
(see |Gun56| lGun76j ). More specifically, given the canonical projection p : H 2 — > T\M 2 
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and the vector bundle ( = (E, T\H 2 , ir, V) over the compact surface, its pullback p*(C) is 
trivial, i.e., the total space is isomorphic to H 2 x V. Now any covering translation 7 G T 
extends to a bundle homomorphism 7 on p*((), and this can be written as 

l(z,v) = (jz } a(z,~/)v) , (3.1) 

where a : H 2 x T -> GL(V). 

For our purposes it will be useful to consider the subgroup f C SL(2,R) that corre- 
sponds to the Fuchsian group T C PSL(2,R) via T = f/{±Id}. Instead of a one hence 
considers a function j : H 2 x T — > GL(V), called factor of automorphy. To define the ac- 
tion of 7 G T unambiguously one therefore requires a suitable character \ of T (multiplier 
system) depending on the value of j(z, —Id). Based on these observations sections in ( can 
be obtained from global sections in p*(()'- any vector- valued function %p : H 2 — > V with 
ipipfz) = x{l) j{ z il) ^( z )i called automorphic form, yields a section in and vice versa. 

Definition 1. Let T C PSL(2, R) be a strictly hyperbolic Fuchsian group, then a multiplier 
system of weight k is a unitary character % of f with x(— Id) = (— l) fc . A factor of 

automorphy for the Maafi-Laplacian of weight k is defined by j y (z, k) := |^^ fc/2 when 

7=(°J). Then ip : H 2 — > C is an automorphic form for the Maafi-Laplacian of weight k, 
if 

4>M = X(l)h(z,k)ip(z) , V7 G f . 

The space of all such automorphic forms is denoted by L #(r, k, x)- Similarly, L £. l (T, k, x) 
and L £ 2 (r, k,x) are the spaces of differentiable and square- integrable automorphic forms, 
respectively. In the latter case, integrability is meant with respect to the measure d/i(z) = 

dx o V over a fundamental domain F C H 2 of T. 

y 

It is well-known that acting on a (twice differentiable) automorphic form of weight 
k yields again an automorphic form of weight k (see |Roe66| IRoe67| |Hej76| ). We denote 
the space of all eigenforms of — with eigenvalue A by l Ja(F, k, x)- in case k > 2 there 
exist certain special eigenvalues that are known explicitly, 

{k) k-2j f k-2j \ . 
j = ^2—\ l -^2—) ' J= ' 1 ' 2 '-' 

where [x] is the integer part of x G R. For multiplier systems x ^ 1 the eigenvalue ()3.2j) 
has a multiplicity (g — l)(k — 2j — 1), see Hej76 ; here g > 2 is the genus of the closed 
surface T\H 2 . 

The relations (|2.1jl and ()2.3j) now suggest the following definition of automorphic forms 
for the weighted Dirac operator. 

Definition 2. Let 

j ( z k \.- (h&k) \ 

w.*;- ^ j 1 (z,k-2)J ■ 



k-i 



(3.2) 
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Then we define the space $(T, k, x) of automorphic forms for the Dirac operator with 
weight k to consist of the functions \I/ : H 2 — > A 2 = C 2 that transform as 

*M = x(i)J 7 (z,k)*(z) , v 7 er. (3.3) 

The spaces € l (T,k,x), £ 2 (T,k,x) are defined analogously, and $ P (T, k,x) denotes the 
eigenspace of — corresponding to the eigenvalue p. 

In order to demonstrate that the weighted Dirac operator maps a differentiable 
automorphic form of the type ()3.3|) to a form of the same type, one has to show that under 
a covering translation 7 the operator behaves as 

D k h+ J 7 (^,A;)D fc J-\z,k) . (3.4) 

This can be done in a straightforward calculation. Hence, the automorphic forms ()3.3|) 
represent sections in spinor bundles over the surface r\H 2 . 

We also note that is essentially self-adjoint on k, x)- The proof is practically 
identical to the case of the Dirac operator (k = 1), which is well-known (see, e.g., |Fri97j 
where also Spin c -connections are considered). Thus we conclude that the spectrum of D^ 
is real and discrete. 

We are now able to establish a connection between the eigenforms of the Laplacian and 
of the Dirac operator. For this we recall that the spectrum of — is bounded from below 
by |(1 - |) |R.oe66| . 

Proposition 1. Let = (l^) be an eigenform of — D& with eigenvalue p, and let X = 
p 2 + |(1 - |). Then ipi G % A (r,fc,x)- On the other hand, if G ^(r, k, x), then * = 
(iA*fi/>) ^ i?p(r, ^, x)- Moreover, if eigenforms of the Dirac operator (p 7^ 0) are linearly 
independent, the same holds for the corresponding eigenforms of the Maafi-Laplacian, and 
vice versa. 

Proof. As both operators are elliptic all eigenforms are smooth. The relation between p and 
A can be obtained by using eq. ()2.3|) . To establish the desired transformation properties one 
has to write down the conditions on the components. The rest then follows from |Roe66j . 
The linear independence can be obtained from some straight-forward manipulations using 

-iA^i + p^ 2 = , 

which can be solved for i/j 2 if P 7^ 0. □ 

4 Symmetries 

Apart from the fact that spectra of weighted Dirac operators on compact surfaces are real 
and discrete, further spectral properties can be concluded from symmetry considerations. 
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Lemma 1. The weighted Dirac operator possesses a chiral symmetry, i.e., if = (^) ^ 
$ p (T,k, x), then (J^ 2 ) £ $- p (T,k,x)- Hence the spectrum o/D& is symmetric about zero. 

Proof. The proof amounts to a simple calculation which can be performed conveniently by 
using the relations in |Roe66j once again. □ 

The symmetry statement extends to zero modes in that the Atiyah-Singer index the- 
orem (see |Hit74t IAS68] ) implies that their number is even. From now on we will denote 
this number by 2N. 

A further symmetry is concerned with quantum mechanical time-reversal. For the Dirac 
operators considered here this is implemented through the anti-linear operator T := io^C 
acting on automorphic forms of arbitrary weight, where C denotes complex conjugation 
and a% = ( 9 q 1 ) is one of the Pauli matrices. Hence T 2 = —Id. Moreover, on the Hilbert 
space £ 2 (r, k, x) this operator is anti- unitary. 

Lemma 2. Let ^ G 5p(r, k, x), then T\l/ e §p(r, 2 — k, x)- In particular, if k = 1 and the 

multiplier system x is real-valued every eigenvalue ofD has a multiplicity of (at least) two. 

Proof. The first part of this lemma can be checked by straight-forward calculations that 
we omit here. For the second part one merely has to use the fact that T is anti-unitary 
with T 2 = —Id. This immediately yields that the eigenforms \I/ and T\& with eigenvalue p 
are linearly independent. □ 

At this point we add a few remarks: 

(i) Above we have always assumed that V contains, apart from unity, only hyperbolic 
elements; hence there are no elements of finite order. Thus one can construct a real- 
valued multiplier system by simply assigning ±1 to the generators of V and extending 
this to r. If one also allows for elliptic elements this clearly will not work. 

(ii) Lemma 121 says that and D 2 _fc possess identical non- vanishing eigenvalues (includ- 
ing multiplicities). According to (|2.2|) the weights k and 2 — k, up to sign correspond 
to the same magnetic field strengths. We will therefore restrict subsequent discussions 
to the case k > 1. 

(iii) In general, quantum systems with half-integer spin share the same behaviour un- 
der time-reversal as in the present case. For the respective quantum Hamiltonians 
the multiplicity two of their eigenvalues due to time-reversal symmetry is known as 
Kramers' degeneracy |Kra30j . Moreover, magnetic fields usually break time-reversal 
invariance. This is the reason for the special role played by the weight k = 1. 

(iv) In his tenfold-way scheme Zirnbauer classified quantum systems according to their 
behaviour under basic symmetry operations like time-reversal, rotations, and chiral 
transformations _Zir96 . The resulting symmetry classes are unambiguously linked 
with Cartan's ten classes of symmetric spaces. In this context the present case, with 
time-reversal symmetry, T 2 = —Id, and chiral symmetry is identified as the type CTI 
and is related to the symmetric space Sp(p + g)/Sp(p) x Sp(g) (of compact type). 
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(v) Following the conjecture of Bohigas, Giannoni, and Schmit [BGS84j one expects that, 
generically, correlations among the eigenvalues of geometric operators on manifolds 
of negative curvature can be described by random matrix theory. (See, however, 
exceptions when the Fuchsian group is arithmetic [BGGS92, BSS92 .) The relevant 
random matrix ensemble for the present case is the chiral Gaussian symplectic en- 
semble (chGSE). 

From supersymmetry considerations (see, e.g., |Tha 92j) it is known that eigenvalues of 
the square of a Dirac operator have even multiplicities. We stress that this observation is 
unrelated to Lemma 121 Rather, this degeneracy stems from chiral symmetry (Lemma 
combined with the squaring. 

However, the degeneracy due to time reversal somewhat unexpectedly extends to the 
spectrum of for odd weight. To this end we recall that the Maafi operator raises the 
weight by two, without changing the eigenvalue of an eigenform of the appropriate Maafi- 
Laplacian |Roe66| lRoe67j . It hence maps ^Sx(^, k, x) to ^xiT, k + 2, x)- The operator 
A k+2 is its formal adjoint from L £ 2 (T, k + 2, x) to L £ 2 (T, k, x)- As a consequence, up to the 
special eigenvalues (|3.2J) the spectrum of — depends only on k mod 2. 

Lemma 3. Let V be a strictly hyperbolic group and x a real-valued multiplier system. If 
the weight k is odd the eigenvalues of the Maafi- Laplacian have a multiplicity of (at 
least) two. 

Proof. For odd k = 2m + 1, the multiplicities of the special eigenvalues ()3.2|) are known to 
be (g — l)2(m — j), see |Hej76 and below equation (j3.2jl . In all other cases the statement 
follows immediately from Proposition ^ and Lemma El □ 

Apart from the special eigenvalues, a twofold degeneracy of the eigenvalues can also 
be obtained in a constructive way. To this end, one first lowers the weight from k = 
2m + 1, m £ No, to one by successive applications of operators A/, then applies complex 
conjugation (i.e., time reversal for spin zero), and finally raises the weight back to k — 
2m + 1. Altogether, one thus applies the anti-linear operator 

L S2m+l '■= CA_2m+lA_2m+3 • • • ^2m+l ■ (4-1) 

Now some straight-forward calculations show that this way an eigenform if) is mapped to 
a linearly independent eigenform <fi = L S2 m +i' l P of the same weight and with the same 
eigenvalue. In particular, <p may only vanish when its eigenvalue is of the special type 

We stress that this argument fails if k is even. In the case k = this is easy to see: 
if if) G L Jp(T, 0, x), then Cip G ^^(r, 0, x) as l° n g as X is rea l- But it is well-known that 
for k = the eigenforms can be chosen to be real-valued, so if) and Cif> are not linearly 
independent. 

Of course, according to Proposition ^ the degeneracy in the spectra of Maafi-Laplacians 
implies a corresponding degeneracy in the spectra of the Dirac operators D^. Except for 
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the eigenvalues that derive from the special eigenvalues ()3.2|) of the associated MaaB- 
Laplacian, one can again devise a constructive approach which reveals more clearly that 
the degeneracy is a consequence of a 'generalised time reversal symmetry'. To this end we 
introduce the operator 

A t _ fp'K k 



ik pK fc _ 2 , 

This maps $ p (T,k,x) into $ p i(T,k + 2,x), where p' is determined by p 2 + k = p' 2 and 
sgn(p') = sgn(p). From |Roe66t Satz 5.3] one can deduce that Aj. indeed is an isomorphism 
between # p (r, k, x) and $ p >(r, k + 2, x) if p 7^ 0. Its formal adjoint reads 



A 



fc+2 



p'Kk+2 ik 
pA k 



The analogue to the operator (j4.1|) now is 

SWn+l := B 2m _i ■ ■ ■ B1TC3 . . . C^m+l , 

where B\ := -A\ and = yA k . A direct calculation shows that S 2m +i does not depend 
on p. 

Proposition 2. Let k = 2m + 1 and assume that the multiplier system x is real-valued. 
Then the eigenvalues of the weighted Dirac operator D 2m +i occur with multiplicities of at 
least two. Moreover, 



>2m+l 



A_2m+3A_2m+5 • • • A 2m +l 

A_2 m +lA_2 m +3 . . . A2m-1 



is identical on all spaces 3p(r, 2m + 1, x) of eigenforms, and can therefore be extended to 
all o/£°°(r,2m + l,x). Then * G $ p (T,2m+l,x) implies <$> := S 2m +i^ e ^ p (V,2m+l,x) 
and ($, \1/)l 2 = 0. For m > the form $ vanishes identically iff p is related to one of the 
special eigenvalues \3. ty) via \j = p 2 + |(1 — |). 

Proof. The statement about the multiplicities of eigenvalues is implied by Proposition ^ 
and Lemma 01 The conclusions $ G ^(T, fc, x) and ($, 1 I / )l2 = follow directly from the 
preceding considerations. One must only take care of the possibility that $ may vanish. 
But evaluating ||$||l 2 shows that this norm is zero, iff A = p 2 + |(1 — |) is of the form 



5 Point-pair invariants 

From the case of the Laplacian it is well known that setting up a Selberg trace formula 
amounts to calculating traces of Hilbert-Schmidt operators L that commute with the Lapla- 
cian jSe!56t |Hej761 |Hej83| . The kernels of such operators can be obtained from suitable 
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Poincare series over so-called point-pair invariants, whose construction will be briefly out- 
lined in this section. 

We first introduce two matrix-valued functions on the upper half-plane, which help to 
study transformation properties under fractional linear transformations, 



A(z', z) 



B(z'z) 



(z'-z)i 




\ 

{z'-z)\ 

( («'-«)* o ^ 
(z'-z)i 

n (z-z')h 



Definition 3. Let $ = (|* : [1, oo) — > M. 2x2 be continuous. Assume moreover that 
$2 = ^3 and that there exists a constant C t such that each component satisfies 



(5.1) 



for some e > 0. Then we call 

K(z', z) = -iA(z', z)$ (±(cosh(d(z', z) + 1)) B{z' , z) , 

where d(z', z) denotes the hyperbolic distance of two points z, z' G H 2 , a point-pair invari- 
ant. 

Applying any transformation 7 G SL(2,R) yields 

K{ 1Z ', 72) = J 7 (z', 1) (y, 2) J" 1 ^, 1) . 

The restrictions on the components of <£> assure that K is hermitian, i.e., K(z',z) = 
K'(z, z'). In order to proceed we need a technical Lemma, which can be found in |Hej 76| . 

Lemma 4. Let f : [1, 00) — > C and e > 0. Moreover, assume that 

holds for all a and a constant C e . Then the sum 

^/(§(cosh(d(z', 7 z)) + l)) 
7er 

converges absolutely and uniformly for all z', zetf. 

Starting from the point-pair invariant we can now construct automorphic kernels via 
Poincare series. Due to our restrictions on $, 



K r (z', z) = \Y, K(z', 1Z ) X {l) Mz, 1) 



(5.2) 



7er 
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is well-defined. Some standard manipulations show that the behaviour of K r under trans- 
formations is given by 

K r ( llZ ', l2 z) =x(li)J ll (z',l)K T (z',z) J~\z,l)x-\l2) ■ (5.3) 

This allows us to define automorphic kernels. 

Definition 4. Let F be a fundamental domain for the Fuchsian group V. We say that 
K r : B 2 x B 2 -> C 2x2 is an element of £ 2 (T\H 2 , T\H 2 , x), if (Q holds for all 7 G f and 

||K r ||L2 := y y tr (i^', «)ifr(z / , z)) djx(z) dp(z') < 00 . 

F F 

It is well-known that such an automorphic kernel defines a Hilbert-Schmidt operator 
L:£ 2 (r,l, X )^£ 2 (r,l, X ) via 



[L*](/) := / K r (^)*(z) dfi(z) 



If i^r is hermitian there exists a basis {^/ n } of ortho normal eigenforms of L in £ 2 (T, 1, x) 
such that K r (z', z) — ^ a n \l/ n (z / )\I/J l (z). Furthermore, if was constructed from a point- 

pair invariant via (|5.2|) . then one can easily check that the operator L has a finite trace 
that is given by a n . In the following section we will construct point-pair invariants from 

n 

a Green's function and then calculate such traces. 



6 The trace formula 

Before we can proceed to introduce point-pair invariants we have to identify Green's func- 
tions for Dirac operators on surfaces r\H 2 . For this, and the following, we restrict our 
attention to the Dirac operator itself, i.e., to the weight k — 1. The starting point will 
be Green's function for D on the hyperbolic plane, from which the corresponding Green's 
function on the compact surface can be obtained in terms of a Poincare series. 

Since the resolvent of D is a bounded operator when Im(p) < 0, we can make the ansatz 

(D + p)- 1 #0') = / G(z', z- p)V(z) dfi(z) , 

with G(-,-;p) : H 2 x H 2 -> C 2x2 . For the matrix entries of G we use the notation G = 
( G3 G4 ) ■ Then G can be uniquely characterised as a solution of the matrix-differential 
equation 

(D + p) G{z', z- p) = , for z^z', (6.1) 

with a specified behaviour in a neighbourhood of z = z'. That is, for the diagonal matrix 
entries 

lim [Gi(z',z;p) - -y- \og(d(z' , z))) < 00, £ = 1,4, 

d(z',z)-*0 \ 47T / 
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is required, whereas the non-diagonal entries are regular in z = z' . In addition, Gi(z', z\ p) 
must approach zero as d(z', z) — > oo. 

Furthermore, the transformation property ()3.4j) implies the corresponding behaviour 

G(z', z; p) = J-\z\ 1) G{ 1Z ', 1Z - p) J 7 (z, 1) 

of the Green's function under an isometry 7 G PSL(2,IR). 

In order to solve for this Green's function we closely follow [Ro e66t [Roe67j . where the 
corresponding problem for the Maafi-Laplacians is treated. We begin with introducing 

G{z',z ]P ) :=A~ l (z',z)G(z',z;p)B- 1 (z',z) , 

which is invariant under PSL(2,R), and note that D := A" 1 (z , i) D A(z , i) is invariant 
under the stability group PSO(2,R) of z = i. Therefore, the differential equation for 
G corresponding to (|6.1|) is transformed into polar coordinates (a, <fi) for z G H 2 , where 
cr = i(cosh(<i(z, z')) + 1). With H(a;p) = G(z',z;p), this leads to a system of ordinary 
linear differential equations for H = ( ^ ^ ) , 



H(a;p) 



p 

U<j(*-l))-hH 3 (a;p) 



Ua{a-l))-^H 2 {a-p) 



The solution of (j6.2j) for a — ► 00 is given by 
Hi (a;p) = 



47T 



r(2ip+ 1) 



1 



--^-«7-2 / ^(1 - tV'-Va - t)- ip dt 

47T 



(6.2) 



(6.3) 



and 



H 2 (a;p) 



~a- l ^{a - 1)I^±11f(1 + ip, 1 + ip; 1 + 2ip; I) 



47T 



r(2ip + 



47T 



(a - 1)3 / _ i)*" 1 ^ _ t)-^- 1 dt 



(6.4) 



where F(a,b;c;z) is a hypergeometric function (see, e.g., |Erd53p . Moreover, H 3 (a; p) = 
H 2 (a] p) and H^a; p) = Hi(a;p). From this representation we infer an upper bound for 
the components of the Green's function, 

\Gi(z', z; p) < const. \p\ ^-MpM*' ,z) ^ if z )>d >Q . (6.5) 
As the singularity of G for d(z', z) — > is integrable we obtain the following Lemma. 
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Lemma 5. Let f : H 2 — > C 2 fre bounded and continuous, then 

G{z',z;p)f{z) dp(z) 



H 2 



converges absolutely and uniformly in Re(p) as /ong as Im(p) < — |. 

With Lemma 0] and the estimate (|6.5|) . Green's function for the Dirac operator on a 
compact surface can now be obtained in terms of a Poincare series (|5.2jl . 

Lemma 6. Let Im(p) < — 1 ; inen 

G r (z / , p) := - GO', 7^; p) x(t) M z > x ) 

converges for z ^ z' mod T, and is the Green's function for D on r\H 2 . 

We are now in a position to construct point-pair invariants as in Definition El through 



00 

9(a) :=i J H(a;p)h(p) dp 



where h is a function as specified below. 
Definition 5. A function h : C — > C that satisfies 

• h is even, i.e. h(p) = h(—p), 

• h is complex analytic in the strip |Im(p)| < (3 for some fixed /3 > | + e, 

• there exists 5 > 0, such that the bound 

|%)| < const. (1 + Re(p))" 2 "' 5 
holds uniformly for all p in the above mentioned strip, 
is called an admissible test function. 

We denote the Fourier transform of an admissible test function by 



OC 



9{U):= 2^ / h ^ e ~ lpU d P- 



All that we have to check that this indeed leads to a point-pair invariant is the condi- 
tion (JHHJ). 
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Lemma 7. Let o — 1 > k > 0, then the components K,i of K are bounded from above by 

|^(2'2)|<C K e- ( ^' 2) , z = l,. ..,4, (6.6) 
with some C K > 0. Furthermore, the limit lim K(z',z) is well-defined. More precisely, 

d(z',2)— >0 

we have 

oo 



tvK(z, z) = — / ph(p) coth(7rp) dp . 

2% J 



Proof. To prove the first part one shifts the integral by — i/3. Using the estimate ()6.5|) on 
Gi and shifting back, the bound ()6.6|) follows immediately. For the second part we note 
that 

oo 

lin \ T~ / P h (p) ^Sl " - !) d P = > 

cr^l+ 47T J 

— oo 

thus we can add this term to Then again shifting the integral by — ie and using MOS66, 
p. 44] one obtains 

oo— ie 

lim Klz', z) = — / p%) [i>(ip) +V(ip+1)] dp , z = l,4, 



d(«',z)-*0 4i7T z 



-oo— ie 



where ^(z) = -4- loglX-z). After substituting p by — p and using the fact that h is even the 
statement follows immediately. □ 

Having showed that K is indeed a point-pair invariant, we can introduce the auto- 
morphic kernel 

K r (z',z) := ^i^', T *) X (7y 7 M) (6.7) 

which by construction is in £ 2 (r\H 2 , T\EI 2 , x). The corresponding Hilbert-Schmidt oper- 
ator on £ 2 (r, 1, x) is called L. 

Lemma 8. Any ^ e 3p(r, 1, x) zs simultaneously an eigenform of the Hilbert-Schmidt 
operator L, 

m(z') = a(pWz') , 

where the eigenvalue A satisfies the equation 

A(p) + A(-p) = 2%) . (6.8) 

Proof. With 

MOO = 2^/ [E j h(p)G(z' n z;p) X (l)J 7 (z,l) dp) dp(z) 
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a standard calculation yields 

oo— i/3 

MOO = ^ / I / GvV,z\pfMz) M*) I V • 

-oo-i/3 \F 

From this we read off the eigenvalue of L as 

oo— i/3 



2vri 7 p' - p 



-oo— i/3 

and this implies (jfi.8j) . □ 

As already noted in Section HJ the spectrum of the Dirac operator on a compact surface 
is real and discrete. Moreover, according to Lemma ^ the spectrum is symmetric with 
respect to zero. We hence denote the non-negative eigenvalues (listed with their respective 
multiplicities) by < po < p\ < . . . , which include N of the 2N zero modes. The Hilbert 
space £ 2 (T, 1, \) therefore has a basis of orthonormal eigenforms ^ n ,±, 

(D± Pn )^n,± = . 

Thus the automorphic kernel possesses a spectral expansion of the form 

oo 

K r { Z ', Z)=J2 ( A (Pn) *„,+ (*')<+(*) + A(-Pn) *„,-(*')<-(*)) > 
n=0 

which immediately yields the spectral side of the desired trace formula. 
Lemma 9. The Hilbert- Schmidt operator L has a finite trace, given by 

tr(L) = 2^/i(p n ) . 

n=0 

Next we need to compute the geometric side of the trace formula, i.e., 

tr(L) = J K r {z,z) dp(z) , 

F 

with the representation (|6.7|) of the automorphic kernel. Following a standard procedure, 
this yields 



tr (L) = ^x(7) £ tr ( /*( 
{7} [ 3 ]ez 7 \r \/ F 



*,7z) J 7 (z,l) d/i(z) , (6.9) 
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where {7} denotes the r-conjugacy classes of 7 G f and Z 7 is the centraliser of 7 in T. 
Next we use the natural pairing of the disjoint conjugacy classes {7} and {—7}, which 
cancels the factor |. Within these pairs we choose the conjugacy classes with tr(7) > 0. 
Moreover, it is well-known that 

D,:= |J 9(F) 

[g]ez-,\r 

is a fundamental set for the centraliser Z 1 C T. Introducing primitive hyperbolic elements 
7 P and their conjugacy classes, one can rewrite ()6.9|) as 



tr(L) = tr / K(z,z) dp(z) 



\iv\ n=l \ Jd i p 



Upon a conjugation with a matrix in SL(2, M.) any hyperbolic element can be brought into 

the Jordan normal form 7 = ( e ^ °; 7 ) • A fundamental domain for the group generated 

by 7 is given by {z e H 2 |l < y < e lj }. Note that J 7 (2, 1) = 1, so all that remains to be 
done is evaluating integrals of the form 

(e l 00 \ 
J J K(z^ n z) dx^\ . (6.10) 

In order to calculate this integral we introduce 

2 nl fx . nl\ 2 

t = cosh h — sinn — , 

2 \y 2 J 

which allows us to perform the integration with respect to y. Inserting the explicit expres- 
sions (|6.3|) and ()6.4|) for the function H, the integral (|6.1UJ1 reduces to 

/(7 n ) = 7^coth^ / ph(p)x 



—00 

00 



J V 2 /^o r(2ip + m + l)m! P 



Interchanging the order of integration and summation and employing the relation 



1 Y 2q j, 2qT(2q + 2m) ( I ^ 



2cosh(fW ^ T(2q + m + l)m! \2cosh(f 

v ^ ' ' 771=0 v ' x v z 
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e -qa 



which follows from 



UR651 1.114] for a > and Re(q) > 0, finally yields 



oo 



'^ = ^hj^^^^ry (6 - n; 



-oo 



Putting together Lemma [7[ Lemma M and equation (|6.11|) then leads to the desired trace 
formula. 

Theorem 1 (Selberg trace formula for the Dirac operator). Let T C PSL(2,M) be 
a strictly hyperbolic Fuchsian group with fundamental domain F of area A(F) and fix a 
multiplier system x of weight one. Moreover, let {p n }^L be the non-negative eigenvalues 
of the Dirac operator D on £ 2 (T, l,x)> including half of the zero-modes. Then, for any 
admissible test function h, see Definition^ the following trace formula holds, 



1? / PKP) C ° thM d ^ + EE^)^d • (6-12) 
47F L ^>n=i 2smh(^) 



m=0 _ OQ { 7p } n=l 

In accordance with Proposition ^ this trace formula is identical to the one for the 
Maafi-Laplacian — Aj on the same surface and with the same multiplier system. We stress, 
however, a difference in the interpretation of the geometric side. Both expressions can be 
viewed as sums over the closed geodesies on r\EI 2 , weighted with the factors x(lp)- F° r 
the Maafi-Laplacian these factors stem from the non-vanishing magnetic fluxes that are 
necessarily present (see, e.g., [CG093 ), whereas for the Dirac operator with weight k — 1 
there is no magnetic field involved. Here the factors xilp) reflect the fact that the classical 
dynamical system associated with the quantum dynamics generated by the Dirac operator 
is not the geodesic flow, but its associated frame flow, see also | J306I . This is analogous 
to Dirac operators on M 3 [BK99J. 

The frame flow (see, e.g., pBG80 ) is a flow in the frame bundle FM over a Rieman- 
nian manifold M consisting of a parallel transport of oriented orthonormal frames along 
geodesies. To be precise, let {ei(p), . . . , e n (p)} be an orthonormal basis of T P M. Then this 
frame is transported along the geodesic determined by e\ with the Levi-Civita connection. 
This flow can be lifted into the spin structure by assigning an element g(p) e Spin(n) to 
every point along the geodesic. This may also be done in the spin representation yielding 
p{g{p)). In the present case, where M = r\H 2 , p(g(z)) must equal p{g{iz)) for every 
7 G f. However, in the induced (trivial) bundle over the hyperbolic plane the transforma- 
tion property 

P{g{iz)) = x{l) J-tiz, 1) p(g(z)) 

applies, see (13.1)) and ()3.3|) . The factors xilp) m equation ()6.12|) therefore reflect the 
spin structure. Hence, although on two-dimensional manifolds the frame flow is not much 
different from the geodesic flow, only the latter is a natural dynamics that can be lifted 
into the spin structure. For n > 3 the situation will be different since in such a case a 
frame flow yields interesting classical dynamics beyond the geodesic flow. 
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7 Some Applications 

A first immediate application of the trace formula, which can be proved in the standard 
way (see, e.g., Hcj76 ), concerns the asymptotic distribution of the eigenvalues. 

Proposition 3 (Weyl's law). Let N(p) be the number of non-negative eigenvalues of D 
which are smaller than p. Then 



p — > oo. 



Another application consists of determining properties of the related Selberg zeta func- 
tion. If Re(s), Re(cr) > 1 the function 



h(p) :- 



p 2 + (s- \f p 2 + (a 



1)2 
2> 



satisfies the criteria of Definition 03 to serve as an admissible test function. In the trace 
formula (jdl2|) it leads to a relation for the trace of a regularised resolvent, 



£ 

m=0 



1 



Pi + (* - |) 2 



+ 
+ 



A(F) 
1 Z'U 



1 



1 



1 Z'h) 



2s-lZ(s) 2o-\Z{p) 
Here Z(s) is Selberg's zeta function which is defined by 



Z(s) := J] II^ 1 - X(7p) e~M*+-)) , Re(s) > 1 . 
{f P } fc=o 

Now proceeding along the lines of |Ste87t I5ar87| we find an analytic continuation of Z(s) 
into the entire complex plane. 

Proposition 4. The Selberg zeta function for the Dirac operator is an entire analytic 
function. Moreover, it can be represented as 

2N 



Z(s) 



(2N)\ 



e ( 2 > 7D e l 2' 2" x 



x 



A(F) OO 
2ir 



n 

m=N 



1 + 



Pi 



(7.1) 



where G is Barnes' double T -function (see War 9 '9f ) and 2N denotes the number of zero 

modes ofD. The trivial zeroes of Z(s) are given by s = —\ — n with multiplicities (2n)~br~ , 
n G N. The non-trivial zeroes are s = ~ ± ip m , with the same multiplicities as the eigen- 
values p m o/D. 
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In ()7.1|) the constant j£> is a generalised Euler constant as introduced in |Ste87j and 
Z( 2N ) denotes a derivative of the zeta function of order 2N. Moreover, due to the Gau8- 
Bonnet theorem, = 2 (g — 1) is a positive integer, where g is the genus of the surface 

r\e 2 . 
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